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Dynamical overlap simulations using HMC 
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We apply the Hybrid Monte Carlo method to the simulation of overlap fermions. We give the fermionic force 
for the molecular dynamics update. We present early results on a small dynamical chiral ensemble. 



1. INTRODUCTION 

It is impossible to run dynamical simulations at 
realistic masses with Wilson fermions, which ex- 
plicitly break chiral symmetry. It is still uncertain 
whether taking the fourth root of the determinant 
in staggered fermion simulations violates locality. 
However, overlap fermions satisfy a lattice chiral 
symmetry, and have no such conceptial difficul- 
ties. The disadvantage with overlap fermions is 
that the calculation of the matrix sign function 
used in the overlap operator requires 0(100) calls 
to the Wilson operator: overlap fermions are com- 
putationally costly! We have developed methods 
which can accelerate the inversion of the over- 
lap operator by more than a factor of 4 |1I2I,'{| . 
bringing simulations on small lattices and at large 
masses to within the capabilities of modern com- 
puters. In this article, we sketch the development 
of a Hybrid Monte Carlo 0] algorithm for dynam- 
ical overlap fermions (see also |fi!7| h 

2. FERMIONIC FORCE 

The massive overlap operator is given by 

D = (1 + /x) + (1 - /i)7 5 sign(7 5 D W /), 

where D\y is the Wilson operator with a nega- 
tive mass parameter, and the bare fermion mass 
is proportional to — /i). The fermionic force 
is the differential of 1/{D^D) with respect to the 
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gauge field (x) . The calculation of the matrix 
sign function has to be accelerated by projecting 
out the smallest Wilson eigenvectors k^') (with 
eigenvalues A') and treating them exactly. The 
eigenvectors can be differentiated using a proce- 
dure analogous to first order quantum mechanics 
perturbation theory. The fermionic force acting 
on a momentum n conjugate to a spinor field <f> 
is F^x) : 
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where e is the sign function and 
fl = (Q-A,)- 1 (l-|^)<V , |). 



(2) 



The delta function in the fermionic force has to 
be treated exactly, because ignoring it will lead 
to large jumps in the hybrid Monte Carlo energy, 
and an unacceptably small acceptance rate. This 
can be done by adding the following terms to the 
leapfrog update: 



U c = e^ n - (AT / 2) C/_(Ar/2), 



n^(Ar/2) = n J _ v (Ar/2)^l 

U+{At/2) = e -iTcTl + (AT/2) Uci 
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where U is the gauge field, r the molecular dy- 
namics artificial time, Uc the gauge field at the 
moment of the crossing (which can be calculated 
to numerical precision using a Newton-Raphson 
procedure), H N the component of the crossing 
normal the A = surface, and the + sign indicates 
the gauge field/momentum after the crossing, i.e. 
after the smallest eigenvalue has changed sign. 
This procedure is area conserving, reversible, and 
conserves energy up to O(Ar). A simple exten- 
sion of this method will conserve energy up to 
0(Ar 2 ) 0. 

We cannot use the momentum update © when 
1 + 4d/(Tr(n^)) 2 < 0. The authors of (i sug- 
gested that in this case, we reflect the momen- 
tum of the potential wall, setting = — II_ 
(a useful analogy is a classical mechanics particle 
approaching a potential wall). We must use this 
"reflection algorithm" for all d < -(Tr(U N )f/i 
(and there will be no change in the topological 
charge), and the "transmission algorithm" out- 
lined above for d > -(Tr(U N )f/i (which will 
lead to a change in the topological charge). 

3. CHIRAL PROJECTION 

It was suggested in a study of the Schwinger 
model [S] that some gain could be achieved by 



projecting into the chiral sector with no zero 
modes (the zero modes have to be accounted 
for by re-weighting when the ensemble averages 
are taken). Because [D^D,^] = 0, we can use 
the operators D^DP± (with P± = |(1 ± 75)) to 
reconstruct the entire non-zero eigenvalue spec- 
trum of the overlap operator. One can show that 
det-D = /i IQI det(D t L»P±), where Q is the topo- 
logical charge, and the sign is chosen so that we 
work in the chiral sector without zero modes. Be- 
cause D^DP± requires only one call to the sign 
function, we can use this to generate a single 
flavour ensemble in half the time we could gen- 
erate a two flavour ensemble. This method has 
two advantages: it will allow more frequent topo- 
logical charge changes, reducing the autocorrela- 
tion length, and secondly there is no exceptionally 
large fermionic force when there is a zero mode. 
However, it will also generate a large number of 
Q 7^ configurations which will be weighted by a 
small factor when we take the ensemble average. 
To preserve detailed balance, we cannot allow the 
gauge field to move into the "wrong" chiral sector 
during the molecular dynamics 0. 

4. NUMERICAL TESTS 

We tested the algorithm described in section 
121 on 4 4 (3 — 5.4 lattices with masses \i = 
0.05,0.1,0.2,0.3,0.4 and 0.5. Figure [T] plots the 
topological charge and dE/dr (E is the HMC en- 
ergy) for each molecular dynamics time step for 
70 trajectories taken from the /j, = 0.05 ensem- 
bles with and without chiral projection. There 
are no large spikes in the energy difference, show- 
ing that our correction step does indeed con- 
serve energy. Secondly, the chiral projected en- 
semble, as expected, has more frequent topolog- 
ical charge changes, and a considerably smaller 
fermionic force when Q ^= 0. 

Tabic Ogives the plaquette, Polykov loop, and 
topological susceptibility for our ensembles. The 
plaquettes and Polykov loops are in good agree- 
ment between the chiral projected and non-chiral 
projected ensembles; the average values of Q 2 
(proportional to the topological susceptibility) 
are not, for reasons which we will outline in a 
later paper [3] . The general behaviour of the pla- 
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Figure 1. Energy conservation plots for the 4 4 
j3 — 5.4 k == 0.225, /i = 0.05 ensembles, with 
(top) and without (bottom) chiral projection. 
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0.19 
0.40 
0.89 
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3.56 



< Q 2 > P 



< S pf > 



0.016(1) 

0.09(3) 

0.16(3) 

0.24(4) 

0.29(4) 

0.32(10) 



0.061(21) 

0.064(18) 

0.048(10) 

0.027(9) 

0.032(9) 

0.005(10) 



0.578(2) 
0.575(2) 
0.568(2) 
0.562(1) 
0.561(1) 
0.554(1) 



0.05 0.19 0.01(1) 0.063(6) 0.576(1) 
0.1 0.40 0.13(3) 0.071(7) 0.574(1) 
0.2 0.89 0.29(4) 0.043(6) 0.569(1) 
0.3 1.52 0.34(5) 0.048(8) 0.562(1) 
0.4 2.37 .40(5) 0.029(9) 0.559(1) 
0.5 3.56 0.45(14) 0.032(13) 0.554(1) 
Table 1 

The plaquette (S p g), Polykov loop (Pi) and en- 
semble average of Q 2 for the 4 4 k — 0.225 en- 
sembles with (top) and without (bottom) chiral 
projection. 



quette and topological susceptibility are in agree- 
ment with our expectations at these large masses 
(our smallest mass is of the same order of magni- 
tude as the strange quark mass) , and the Polykov 
loop is small, suggesting that our configurations 
are confined. 

5. CONCLUSIONS 

We have shown that it is possible to run dy- 
namical simulations using overlap fermions, and 
that the delta function in the fcrmionic force can 
be overcome. The results on our small test con- 
figurations are sensible. We hope to begin sim- 
ulations on 12 and similar lattices in the next 
year. 
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